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The finite size behavior of the susceptibility, Binder cumulant and some even moments of the 
magnetization of a fully finite 0(n) cubic system of size L are analyzed and the corresponding 
scaling functions are derived within a field-theoretic e-expansion scheme under periodic boundary 
conditions. We suppose a van der Waals type long-range interaction falling apart with the distance 
r as r~ ( - d+a \ where 2 < a < 4, which does not change the short-range critical exponents of the 
system. Despite that the system belongs to the short-range universality class it is shown that above 
the bulk critical temperature T c the finite-size corrections decay in a power-in-L, and not in an 
exponential-in-L law, which is normally believed to be a characteristic feature for such systems. 



I. INTRODUCTION 



It is well known that the critical properties of a given bulk system depend on a small number of parameters like its 
dimensionality, the symmetry of the order parameter and the long-rangeness of the interaction in the system under 
consideration. If the Fourier transform of the interaction v(q) has a small \q\ expansion of the form 

v{q)=v + v 2 q 2 +v a q a +w{q), (1.1) 

with w(q)/q a — > when q — > and a > 2, then the thermodynamic critical behavior of the system is supposed to be 
like that of an entirely short-ranged system |5J]. In the opposite case, when a < 2 the critical behavior differs essentially 
|l],^| from that of the short-range system and is characterized by critical exponents that do depend on a (below the 
corresponding upper critical dimension that is d u = 2a in this case) . On the basis of the above bulk picture one 
normally supposes that in the finite systems the same general property will take place: if a > 2 the finite-size behavior 
will be that of the corresponding short-ranged finite-size systems || , characterized by exponentially fast decay of the 
finite-size dependence of the thermodynamic quantities (at least when the critical region of the system is leaved in the 
direction towards higher temperatures; the low-temperature behavior depends on additional features like existence, 
or not, of a spin-wave excitations - Goldstone bosons). It turns out that the last is not true, at least for 2 < a < 4, 
and an evidence about that within the framework of the mean spherical model has been reported in 0]. For example, 
it has been demonstrated that the finite-size dependence of the susceptibility in such a system is given by (2 < d < 4, 
2<cr<4, d + o- <6) 

X (t, h- L) = L-y/»Y( Xl ,x 2 , bL 2 —^), (1.2) 

or, equivalently, 

X (t, h; L) = W> v [Y sr (x u x 2 ) + bL 2 —W lr (x 1 ,x 2 )] , (1.3) 

where x\ = citL 1 ^ , x 2 = C2hL A / u , and Y, Y sr and Y lr are universal functions (recall that 77 = for the short-range 
spherical model). The quantities ci, c 2 and b are nonuniversal constants. In the high-temperature, unordered phase, 
where tL 1 ^ — > 00, one observes Q that the long-range portion of the interaction between spin degrees of freedom 
gives rise to contributions of the order of bL~( d+cr \ In other words the subleading long-range part of the interaction 
gives rise to a dominant finite-size dependence in this regime which is governed by a power-in-L law. More explicitly, 
one obtains Y sr (xi,0) ~ exp(— const, x"), while 

Y lr ( Xl ,0)~ x - d »- 2 -<, (1.4) 
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when x\ — » oo B. This asymptotic is supported fr om t he existing both exact and perturbative results for models 



with leading long-range interaction |5|-|8[. Note that (1.4) implies for the temperature dependence of this corrections 
that 

X (t, h; L) - X (t, h; oo) ~ t -dv-2 7 L -(d+a) ^ tL i/u ^ ^ ( L5 ) 

In addition, let us note that the standard finite-size scaling [p[-p^[ is usually formulated in terms of only one reference 
length, namely the bulk correlation length £. The main statements of the theory are that: 

i) The only relevant variable in terms of which the properties of the finite system depend in the neighbourhood of 
the bulk critical temperature T c is L/£. 

ii) The rounding of the phase transition in a given finite system sets in when L/£ = 0(1). 

The tacit assumption is that all other reference lengths will lead only to corrections towards the above picture. As 
it is clear from Eqs. ( |1.2| ) - ( |l.5| ) this is not the case in systems with subleading long-range interactions. This is an 
important class of systems. It contains all nonpolar fluids where the dominant interaction is supposed to be of van 
der Waals type, i.e. of the type given by Eq. ( |l.l[) with d = a = 3. 

In fact a similar problem has been recently studied by Chen and Dohm Jlij-p^l . They considered a field-theoretical 
model with short-range interactions and wavelength-dependent cutoff of fluctuations A. They observe corrections to 
the infinite system thermodynamic behavior going as an inverse power law in L that do depend also on LA and not 
only on L/£. As it has been clarified in Q] the power law contributions to the finite size corrections result there from 
the interplay of two features of that model. The first is a sha rp c utoff of fluctuations in momentum space and the 
second is the removal of all the terms beyond the q 2 one in ( |l . l| ) , which has the effect of introducing an effective 
interaction that falls off as a power law in the separation between degrees of freedom. This power-law interaction 
leads immediately to power-law contributions to the finite size corrections. 

Theoretically the critical properties of finite-size systems have been studied on the examples of exactly solvable 
models, by renormalization group calculations - both in the field-theoretical framework and in the real space, by 
conformal invariance and by numerical (mainly Monte Carlo) simulations. An essential part of these investigations is 
well described in a series of reviews pj]- |r^ , [i7| , |T8[ . 

The 0(n) models are the most often used examples on the basis of which one studies the scaling properties of finite- 
size systems. The best investigated cases are those of the n = 1 (Ising model) and the limit n = oo, which includes 
the spherical model JlJ,[l3| ■ The last model is especially suitable for the investigation of its finite-size properties since 
it is exactly solvable for any d even in the presence of an external magnetic field H. For n ^= 1, oo there are no exact 
results and the preferable analytical method for the derivation of the properties of the corresponding models (like 
XY , i.e. n = 2, and Heisenberg, i.e. n = 3) is that one of the renormalization group theory. An important amount 
of information for such systems is in addition derived by numerical simulations, normally via Monte Carlo methods. 
As a rule the investigations are concentrated on interactions of finite range. As examples of long-range interactions in 
addition to the equivalent neighbors the case of power-law decaying interactions have been considered. In the case of 
a < 2 analytically only the finite-size scaling properties of the n = oo limit are well established. For finite n a limited 
number of recent numerical results |l^-|2^], as well as few theoretical works 0,||,|22|,^3| are available. In |23| one studies 
the crossover from long to short-range interactions, i.e. the limit a — > 2 and concludes that the renormalized values 
of the temperature and the coupling constant are continuous functions of a. The case of a < 2 has been inve stig ated 
in references (under periodic boundary conditions). It has been found that, as for the bulk systems flilpif, the 

critical behavior depends on the small parameter e = 2a — d, where 2a corresponds to the upper critical dimension 
in such systems ppif . The results are obtained in powers of \fe. The quantities of interest have been the shift of 
the critical coupling, the susceptibility and the Binder cumulant B at the critical temperature T c |^2],[7j and above 
it as a function of s. It has been found that the numerical results obtained in |22| for the Ising model do not 
agree with the predicted (up to one loop order) behavior of B j7],^2| . In |2lJ one even reports disagreement with the 
well established theoretically fact that the critical exponent of the system do not depend on a if a > 2. In (T^| a 
Heisenberg model and in |2(j an Ising model have been investigated by Monte Carlo methods in the case when their 
critical behavior is characterized by classical (mean-field) critical exponents. 

In the present article we will consider the case of long-range power-law decaying interaction characterised by a > 2 
in its Fourier transform. As it was already mentioned above the recently obtained results for n — oo limit indicates 
that the well-spread opinion that such an interaction is uninteresting for the critical behavior of the finite system || is 
not fully correct. Here, following the method used in JtJ we will generalize the results available for n — oo to the case 
of finite n. We will use e-expansion technique up to one loop order in the interaction coupling. We will investigate 
the behavior of the Binder cumulant, susceptibility, and some more general even moments of the order parameter. 

The plan of the article is as follows. In Section [IH we review, briefly, the <^ 4 -model with long-range interaction and 



discuss its bulk critical behavior. Section III is devoted to the explanation of the methods used here to achieve our 



analysis. We end the section with the computation of some thermodynamic quantities of interest. In Section IV we 
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discuss our results briefly. In the remainder of the paper we present details of the calculations of some formula used 
throughout the paper. 



II. GENERAL CONSIDERATIONS 



In the vicinity of its critical point the Heisenberg model, with short as well as long-range interaction decaying in a 
power-law, is equivalent to the d-dimensional C(n)-symmetric model 



PH{<P) = IJ d d x (V^) 2 + 6(y/V 



2,1 4 

r tp + -u ip 



(2.1) 



where tp is a short-hand notation for the space dependent n-component field ip(x), tq = r$ c + to (to oc T — T c ) and 
Mo are model constants. V is the volume of the system. In equation ( |2.l| ), we assumed % = fcs = 1 and the size scale 
is measured in units in which the velocity of excitations c = 1. We note that the second term in the model denotes 
qr CT |(/3(g)| in the momentum representation where the parameter a > ( with a/2 being noninteger) takes into account 
the contribution of the long-range interactions in the system. In (2.1) j3 is the inverse temperature. Here we will 
consider periodic boundary conditions. This means 



{x) = L d ^2ip(q)exp(iq-x), 



(2.2) 



where q is a discrete vector with components qi — 2irn,i/ L (rij = 0, ±1, ±2, • • • , i = 1, • • • , d) and a cutoff A ~ a -1 (a 
is the lattice spacing). In this paper we are interested in the continuum limit, i.e. a — * 0. As long as the system is 
finite we have to tak e int o account the following assumptions L/a —* oo, £ — > oo while £/L is finite. 

The Hamiltonian ( [2.l| ) is, of course, well known in the literature. First, it has been used to investigate the critical 
behaviour of systems with reduced space dimensionality exhibiting phase transitions ]lj] . Let us recall that in such 
systems a phase transition can occur only if the interaction is long-ranged enough. The critical behaviour of the model 
depends strongly upon the nature of the interaction controlled by the parameter a. With a < 2 it has been used for 
detailed investigation of the critical behavior of 0(n) models including questions like the a, d and n dependence of 
the critical exponents and critical amplitude ratios, as well as for calculation of their values, and for determining of 
the universal scaling functions of both the infinite, as well as of finite systems. In this case the critical exponents of 
the system are a dependent. By increasing a, a crossover from long-range critical behavior to short-range one takes 
place. The crossover happens at a point, which can be determined from general considerations (see for example page 
71 of reference p5|). This 'critical' value of a is given by a = 2 — ry, where rj is the Fisher exponent for the short 

range model. When a > 2 one usually considers the model as equivalent to a = 2 case and omits the b (V 7 / 2 ^) term 
in the Hamiltonian, since it was believed that this term does not contribute to the critical behavior of the system. 
Indeed, in this case, the critical exponents does not depend on the parameter a. As it was already mentioned, such a 
procedure can lead to incorrect results for finite-size systems. This was demonstrated in [|| on the example of n = oo 
model. In the current article we will demonstrate that the same remains true also for a finite n. 

for the case a > 2 is achieved by considering the 



The investigation of the critical phenomena of the model ( 2.1 



long range interaction as a perturbation to the short range one 1 2q-|28| . This allows the adaptation of the theory of 
Feynman diagrams to systems with subleading long range interaction. As a consequence the upper critical dimension 
remains unchanged by that interaction and the critical exponents are those of the model with pure short range 
interaction. The interested reader can find more details in refrences [Efl-GSl . 



Before starting to explore the scaling properties of the field theoretical model (2.1) confined to a finite geometry and 
under periodic boundary conditions, we will give a brief heuristic derivation of the finite scaling hypothesis, based on 
the idea of renormalization group. Here we are interested in the continuum limit when the lattice spacing completely 
disappears. Using dimensional regularization the integrations over wave vectors of the fluctuations are convergent 
and are evaluated without cutoff. When some dimensions of the system are finite the integrals over the corresponding 
momenta are transformed into sums. Since the lattice spacing is taken to be zero, the limits of the sums still extend 
to infinity. 

^From general renormalization group considerations a multiplicatively renormalizable observable X, the suscepti- 
bility for example, will scale like 



X\b,g,b,n,L\ = C,{p)X [t(p),g(p),b(p), /j.p,L] 



(2.3) 



where t = (T — T c )/T c is the reduced temperature, g is a dimensionless coupling constant and L is the finite-size scale. 
The length scale p is introduced in order to control the renormalization procedure. Here b(p) is an irrelevant from RG 
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point of view variable whic h mimics the influence of the subleading long range interaction on the critical behavior of 
the system. Equation ( [2.3|) is obtained using the assumption that the size L of the system does not renormalize [ ^9| . 

It is known (see, e.g., [29]) that in the bulk limit, when g(p) approaches its stable short-range fixed point g* of the 
theory, then we have 



t(p) Ktp 1 '"- 2 , C{p)~P~ u/v ~ P:l and b(p) « bp"*"- 2 , 



(2.4) 



where j x and v are the bulk critical exponents measuring the divergence of the observable X and the correlation length, 
respectively, in the vicinity of the criti cal p oint and p is a scaling parameter. The exponent p x is the dimension of 
the observable X, defined in equation (2.3). The critical exponent 7^ = 2 — 77 — a |3(J. Using dimensional analysis 
together with equation (2.2) one gets 



X[t,g,b,fi,L] = (wy*({p)X [t(p)(pp) 2 ,g(p),b(p)(pp) 2 ,l,L/pp] 



(2.5) 



Choosing the arbitrary parameter p = we obtain our final result for the scaling form of an observable X in the 
case, when there are subleading long-range interaction in the finite system 



X[t, g, b, p, L] = lW"/ (tL 1 /", bL 2 -°-» 



(2.6) 



Here the function f(x) is a universal function of its argument. Note that equation (2J3) is the an alog , for finite system, 
of the result obtained in |3C| ]. In the remainder of this paper we will verify the scaling relation (2.6) in the framework 
of model (O). 



III. FINITE-SIZE ANALYSIS 



The method we will adopt here is widely used in the exploration of the scaling properties of finite systems in the 
vicinity of their critical point. It is based on the idea of using a mode expansion, i.e. one treats the zero mode of the 
order parameter, which is equivalent to the magnetization, separately from the higher modes (<j ^ 0). The nonzero 
modes arc treated perturbatively in combination with the loop expansion. The finite modes are traced over to yield 
an effective Hamiltonian for the zero mode: 

exp [-Ties] = Tr^ exp [-Ti {<j) q=0 , <j> q ^o)] ■ (3.1) 

After performing this operation one ends up with an effective Hamiltonian of the form (see Appendix |a]) 

H cS = l -L d (r<$> 2 + hj<$>^ , (3.2) 

where the effective coupling constants are given by 

R = r Q + (n + 2)u L- d ^2 ,J, U i CT > ( 3 - 3a ) 
^ r + q 2 + b\q\° 

U = u -(n + 8)u 2 L- d , - * ^2 - (3.3b) 
^ (r + q 2 + b\q\°) 

In the remainder of this paper we will compute, to the lowest order in e = 4 — d, the effective coupling constants, 
with the initial coupling constants renormalized as in their bulk critical theory, since it has been shown that to the 



one loop order the renormalization of the finite theory is a consequence of the renormalization of the bulk one 1 29 



Simple dimensional analysis shows that the effective coupling constants should have the following scaling forms: 

R = L q - 2 fn. (tL l ' v , bL 2 - a -A and U = L d ~ i+2ri fu (tL 1 ^, bL 2 - a - r ') , (3.4) 



for t > 0, where /r and fjj are scaling functions which are properties of the bulk critical point. They are analytic at 
t = 0. This is a consequence of the fact that only finite modes have been integrated out. 

After evaluating the explicit forms of the functions fa and /[/, we can deduce results for the different thermodynamic 
quantities and the expressions of their respective scaling functions. 
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In order to investigate the large scale physics of the finite system, one has to calculate thermal averages with respect 
to the new effective Hamiltonian defined in (3.2). They are related to the thermodynamic functions of the system 
under consideration. The averages of the field <j) are defined by 



mP\_ fd n <f> 4> 2p exp(-H cti ) 



J d"0exp(-H of f) 



With the aid of the appropriate rescaling $ = (UL d ^ 1 ^ 4 <fi, we can transform the effective Hamiltonian into 

H cS = l -z$> 2 + i$ 4 , (3.6) 

where the 'scaling variable' z = RL d / 2 U~ 1 ^ 2 is an important quantity which has been used in many occasions in 
the investigations of finite-size scaling in critical systems (see for example references [ ^lp^ ]). Explicit expressions for 
some thermodynamic averages of the type (3.5) as well as their asymptotics are presented in Appendix [B|. 
With the effective Hamiltonian (|3.6|), we obtain the general scaling relation 



= L~^- 2 ^ L up/2 f 2p [rL 2 - ^ ^ ) 



for the momenta of the field 4>. Having in mind Eqs. ( [3.4| ), we can write down Eq. (3.7) in the following scaling form 

M 2p = L-' p{d - 2+ ^T2 P {tL 1 l\bL 2 " y - ri ), (3.8) 

in agreement with the finite-size scaling predictions of ([2.6|). In equation ( |3.8| ), the functions Ti P (x) are universal. 

All the measurable thermodynamic quantities can be obtained from the momenta M.2 P - For example the suscepti- 
bility is obtained from 

x = - [ d d x (<p(x)<p(0)) = L^^itL^^bL 2 -^). (3.9) 
n Jv 

Another quantity of importance for numerical analysis of the finite-size scaling theory is the Binder's cumulant defined 
by 



1M 



B = l--^. (3.10) 

In the rema inde r of this section we concentrate on the computation of the coupling constants R and U of the effective 
Hamiltonian ( |3.2| ) for the system with subleading long-range interaction decaying with the distance as a power law. 
As a consequence we will deduce results for the characteristic variable z = RU~ 1 / 2 L 2 ~ r, ~ e / 2 , the susceptibility % and 
the Binder's cumulant B. 



A. Computation of the effective coupling constants 

The finite-size corrections to the coupling constant r$ in the mode expansion reads 

R = r + (n + 2)u L- d V 1 (3.11) 

One of the delicate problems in the finite size-scaling theory is the analysis of the sums appearing in the mathematical 
equations, which forms the basis of the investigation of the scaling as well as thermodynamic properties of the system 
under consideration. In our case this means that we have to find a way to evaluate the sum appearing in the right 



hand side of (3.11). In the absence of the long-range interaction term (b = 0) several methods have been developed 
in order to investigate this sum. When 6^0, i.e. in the presence of the non analytic term in q, a step towards the 
solution of this problem has been made in reference Q. It is based upon the idea that in the long distance physics 
one retains only those contribution to the behavior of the quantities involved that are associated with the effects of 
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long-range fluctuations. In other words we will consider the leading behavior that is due to the small q contributions. 
Expanding in q>, we obtain 

R = r Q + (n + 2)u S L (d, r , 2) - (n + 2)u b (l + r Q -£A S L (d, r , tr), (3.12) 

where 

s ^^)=^!S- (3 - i3) 

q^O y 

In order to evaluate the finite-size corrections to the bulk system we have to analyze the finite-size behavior of the 
function Sl{(1, r, a). This is achieved by making use of the identity 



q 2 P 



exp[-{q 2 + r)t]r p 1 *(~p,~rt)dt, p < 1, (3.14) 



r + q 2 j Q 

where 7* (a, x) is a single- valued analytic function of a and x, possessing no finite singularities J33| 

±^r(a + n + l) r(o) ±^L (a + n)n\ 



Identity ( [Dg ) can be proven by integrating by parts the corresponding series representations of 7*. Similar identity 



has been used in J23| for the investigation of the finite-size behavior of 0(n) model system with a crossover from leading 
long-range interaction to the short-range case, i.e. a — > 2~ 
With the help of this identity one obtains (0 < p < 1) 



1 f n 2 P 

S L (d,r,2(p + l)) = ^ J dq -3_ + L^X^ s irL\dl 



where 



with 



■nyl/2 



(3.16) 



du, (3.17) 



OO j 



Finally for the effective coupling constant R we obtain 

R = ro + (/ ~ b I + ( " + 2 ^ L2 ~ dl -^ ^ L ^ d ) 

-(n + 2)u<,& + r A) L 4 ^-/^ (r L 2 , d) , (3.18) 

Now, we renormalize the theory by introducing the field theoretical renormalization constants, i.e. the scale field 
amplitude Z, the coupling constant renormalization Z g , and Z t - renormalizing the ip 2 insertions in the critical theory. 
This allows to replace the model bare constants r$ and uq in the last equation by their renormalized counterparts 
trough: 

t = ZZ^ 1 (r -r 0c ) and g = ^ Z 2 Z~ 1 u Sj 1 , (3.19) 
where /1 is a renormalization scale, which will be set equal to 1, Sd — |(4'7r) d ' 2 r(d/2) is a phase space factor and 



Z = l + 0(g 2 ), (3.20a) 
77 + 2 

Z t = l + —— g + 0(g 2 ), (3.20b) 

Z g = l + ^g + 0(g 2 ), (3.20c) 



G 



are the usual renormalization amplitudes to one-loop order. 

Finally, using dimensional regularization, at the fixed point g* 
we obtain 



BIS = y[l 



-eS 4 



en + 2 



2n 



In V 



e n - 



4n + 8 



! &L 2- CT ( 2 + cr ) 7r „<t/2 



+ 0(e 2 ) of the theory, in the case d + a < 6 
n + 2 



-I) 



■ t - 



n + 8 



&4 ^scaling (V ') 4) 



^L 2 



1 + 2/ 



scaling 



(3.21) 



2 - ^p§£ + O (e 2 ) . Equation ( |3.21| ) shows that 



where we have introduced the scaling variable y = tL^I" with v 
the effective coupling constant R has the scaling form predicted in equation (3.4). At this order, the exponent n = 0, 
and verifying t he po wers of 77 in the this expression requires a higher order computation. In the particular case 6 = 
from equation ( 3.21 ) we recover the result of reference [pl| . 

When the system under consideration is confined to a finite geometry, instead of the coupling constant uq, we have 
the shifted effective coupling constant U given by: 



U = u - (n + 8)ulL- d 



q#0 



(ro 



b\ q \°y 



(3.22) 



Remark that the summand in the right hand side can be expressed as the derivative of the summand in right hand 
side of Eq. (3.11) with respect to ro. Consequently the result for the effective coupling constant U can be derived 
easily from that of R. Using that observation one gets 



UL e 



-5 4 (l + |(l + rny)) + 



r 2 bL 2 - a „ aw(a + 2) , 



f 2 8 
'n + 8 4 dy scaUng ' 



n + 8 



S\bL 2 - G 



8 



d 2 \ r z=* 



' dy dy 2 



scalins 



(y, 4), 



(3.23) 



at the fixed point, in agreement with the scaling relations of equation (3.4). Equation ( p.23j ) generalizes the resuls of 
refrence ]3l| ] to the case when subleading long range interaction is taken into account. 

Note that the effective coupling constant U has a finite limit at the critical point, i.e. in the limit t — * 0. Indeed, 
as the reduced temperature vanishes it is possible to use the expansion 



-^scaling (V,4) 



scalii 



(0,A) + ^y(C-lny)+O(y 2 ), 



(3.24) 



where 



exp 



u \ u 2 , A . u 2 

■77 - —xA (u) + — 



7T- 



2.2064. 



(3.25) 



After substitution of (3.24) in (3.23) the terms proportional to logy cancel, which shows that the coupling constant 
U is finite at t = 0. Whence, one gets (for y — > 0) 



UL £ = 



n + 8 



:Sa 1 



c - 



9 a 2 
> 1 

<9y 9y 2 



scalint 



(y.4) 



(3.26) 



y=0 



showing that it is possible to evaluate U at the critical point, i.e. it is safe now to set y = 0. 



B. Some thermodynamic quantities 



1. Binder's cumulant 



In this subsection we are interested in the evaluation of the Binder's cumulant ratio, which plays a fundamental role 
in the investigation of the finite-size scaling theory by numerical means. Here we will give only the analytical expres- 
sions. Unfortunately there are no numerical simulation which can approve or not the results we obtain throughout 
this paper. 



7 



Close to the critical point, i.e. in the region tL 1 ^ < 1, we obtain for for the Binder's cumulant ratio 

r 2 [H L 



8 = 1- 



12I*[£(n + 2)] \ 

/ r[l(n + 6)]r[|(n + 2)] r 2 [l(n + 4)] 
I r[I(n + 4)]r[|n] + V[i(n + 2)] 



r[|(n + 6)] r[l(n + 2)] _ r [|(n + 4) 
r[I(n + 4)] + r[Jn] " 




(3.27) 



The cumulant 8 is a function of the variable z, which is itself a function of the scaling variable y. So, a knowledge 
of a final expression for the function z, which appears in the all thermodynamic quatities, is enough to deduce the 
value of the Binder's Cumulant. At the fixed point, we obtain (for y « 1) 



RL 2 



VTn? 



n - 



eS 4 
n + 2 



fixcdpoint 

y 



■ in y 



in e (2 + tr)7r a/2 



n + 8 
1 



eS 4 



-^scaling 



6L 



2-cr 



16 sin 1 7r- 

V ( 



dy 



8 



2-<T 



1 + 2/ 

» | . 

9y 3y 2 



scalint 



(y, 4) 



a 2 \ ^2 



scalint 



(2/,4) 



(3.28) 



This expression shows that the Binder's Cumulant B has the required scaling form. At the critical point, i.e. at y = 0, 
we get 



**(0) 



-4V2 



7T y/n + 

This result is obtained with the help of the formula 

2 



In 2 + 6L 2 -'(27r)'- a (l - 4*- x )C (l - \ ) C (2 - |) 



duu 



A\u) - 1 - - 
v it 



= 8(1 - 4 1 ^) 7 r 2 ( 1 ^)r( I ,)C(^ - 1)CM, «/ ^ 0, 2. 



(3.29) 



(3.30) 



Equation ( 3.29| ) is a generalization of the result of j3lj obtained for the model with pure short range forces. Note 
that the form of the expansion in terms y/e is kept but the coefficient is altered and now it is a function the parameter 
<7 controlling the long-range interaction. 

Now we turn our attention to the behavior of Binder's cumulant ratio in the limit z» 1. In this case we obtain 
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--) 


H +o (?)] 




n) 





(3.31) 



wherefrom one has B n {oo) = |(1 — 1/n). This result corresponds to a n-dimensional Gaussian distribution for 
n independent components $].,••• , of the vector variable <I>. For such a distribution it is easy to show that 
M 2 = n < <& 2 >, and M4 = n < $f > +n(n — 1) < $ 2 > 2 , where <£>i is any of the components of the vector 
and < ■ • • > means average with respect to one-component Gaussian distribution Q\. Having in mind that for Q\ 
< $| >= 3 < $ 2 > 2 , one directly obtains that B n = |(1 — 1/n), in a full agreement with the above renormalization 
group result. Obviously, all limiting values lie in the interval from 8 = (Ising model, n — 1) to 8 = 2/3 (spherical 
model, n = 00). 



2. Magnetic susceptibility 

The system we consider here is confined to a fully finite geometry. In this case it cannot exhibit a true phase 
transition, i.e. the thermodynamic functions are not singular. In the vicinity of the critical temperature, which 
corresponds to the region y -C 1, the susceptibility behaves like 
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x = - 




r|j(n+2) 

r[f] 



(3.32) 



The susceptibility in this case is analytic at i = 0. This is a consequence of the analyticity of the effective coupling 
constants R and U. In order to get the final result for the susceptibility one has to replace R and U by their respective 
expressions. After performing this we find that x has an expansion in powers of y/e. This results is valid as long as we 
are concerned by the case d + a < 6. Once we have d + cr = 6,alnL will appear in the expression of the susceptibility. 



The source of this InL is coming from equation (C3) for the coupling constant U at the critical point (see Appendix 
^). This is an extension to finite n, by means of perturbation method, of the result obtained in reference [Q for the 
spherical model. 

In the region corresponding to y 3> 1, we have 



X 



1 

R 



1 - 



O 



1 



(3.33) 



Substituting the effective coupling constants R and U by their respective expressions from equations ( 3.21 ) and ( 3.23 ) 
and using the asymptotic expansions of Jf ca i; ng derived in [Q, we get 



X = Xo 



2S 4 



bL 2 -° f C 4 ^j/- 2 



,*/2 



4S 4 sin7T(r/2 



where 



(l+p)4 1+ fr(I+p + d/2) 



T d/2 



n-p) 



fc#0 



1 



(3.34) 



(3.35) 



Expression ( 3.34 ) for the susceptibility shows that it has the form given by the scaling hypothesis (3.9). It demonstrates 
also that in this regime the critical properties of the system are dominated by the bulk critical behavior, with finite-size 
corrections in powers of L. 



IV. DISCUSSION 



In the present article we have investigated the finite-size scaling behavior of a fully finite 0(n) system with periodic 
boundary conditions and in the presence of a long-range interaction that does not alter the short-range exponents of 
critical its critical behavior. The small \q\ expansion of the Fourier transform of the interaction v(q) is supposed to 
be of the form 



v{q) =v a + v 2 q 2 + v a q a + w(q), 



(4.1) 



with w(q) — > 0/q a , when q — > and 2 < a < 4. In the real d-dimensional space one can think about interactions 
decaying as r ~( d+a K This is an important class of interactions that include also van der Waals type interactions. 

For such a system, in the present article we have demonstrated that all the even moments of the magnetization 
M.2pi including the susceptibility, can be written in the form 



M2 P = L-^ d - 2+r >^ 2p (tL 1 ^,bL 



2-<r-m 



(4.2) 



(see Eqs. @, ( |3~2l| ), ( p3; ), ( |3~32l) , (|3~3^)). Note that one has two scaling variables needed in order to describe in a 
proper way the finite size behavior of these quantities. A special attention has been paid to two important quantities: 
the Binder's Cumulant and the susceptibility. 

In the region away from the critical point we obtained for the Binder's cumulant ratio the expression 



B 



1 / 2 
1-r 1+- 
3 V n 



(4.3) 







with finite size correction falling off in a power law. The above result corresponds to a n-dimensional Gaussian 

distribution for n independent components of the vector variable. Obviously, all the values lie in the interval from 
B = (Ising model, n — 1) to B — 2/3 (spherical model, n = oo) . 
For the susceptibility, when tL 1 '" 1, one has (see Eq. (3.34)) 



X = Xc 



1-e 



71 + 25*4 
n + 8 y 



bL 



2-a 



C^ a^iy 2 



W 2 



45*4 sin7rcr/2 



(4.4) 



One observes that in this regime the susceptibility approaches its bulk value not in an exponential-in-L, as it is 
usually believed to be the case for systems with short-range critical exponents, but in a power-in-L way. The last 
goes beyond the standard formulation of the finite-size scaling, but is completely consistent with the intrinsic large- 
distance power-law behavior of the correlations in systems with long-range interactions (see, e.g. |35[ and references 
cited therein). 

Since r\ = 0(e 2 ) in 0{n) short-range models, we were unable to verify the predicted dependence of the scaling 
functions on 77, which requires calculations up to second order of e, while we have retained only corrections up to the 
first order in e. We hope to return to this problem in the future. 
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APPENDIX A: CONSTRUCTION OF THE EFFECTIVE HAMILTONIAN 



Let us start with the bare Hamiltonian ([2.1 



H{<p} = - J^d d x 



(Al) 



where the spatial integration is over a system of linear extent L in each of its d dimensions. The partition function is 
given by 



Z = J T>(pexp(-H). 



(A2) 



Following reference R2] , we spilt the field 



ip{x) = (f> + X 



(A3) 



into a mode independent part <p, which defines the magnetization, and a part depending on the nonzero modes 
S = L~ d X) q ^o vio) ex P (iQ ' x )- For further calculation we introduce the auxiliary Hamiltonian 



(A4) 



and we treat the rest of the Hamiltonian by using perturbation theory. Within this approximation the partition 
function reads 



Z = 



J 2tyexp(-iToM-r(0), 



(A5) 



where 



r=-ln / X>£cxp(-H(0,£)+W o (<£))- 



(A6) 



Writing the difference between the bare Hamiltonian (Al) and the auxiliary Hamiltonian (A4) in the form 
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H{<j),Z}-H {<jy} =i J d d x[{r + 3u (t> 2 )T, 2 + (V£) 2 + &(V /2 s) 2 ] (A7a) 
+ i,io^d d a;[20I] 3 + il] 4 ]. (A7b) 

keeping in mind that the additional term involving J v d xYl vanishes one gets, after some straightforward calculations, 
including the evaluation of the integrals over the field S, 



T(0 2 ) = ln t r ° + q2 + bq ^ + V n + ^o^i^iCro, L) - i(n + 8) U 2 4 L d 5 2 (r o , L) + • • • , (A8) 

q#0 

where 

5 m (r ,L) =L- d Y - } - r , (A9) 

and the dots represent terms with higher order in (j>. 

Substituting expression ( |A9| ) into that of the the partition function (A5), we end up with the final expression for 
the effective Hamiltonian 

H cS . = ^L d [Rtf + hj<j>% (A10) 
where the effective coupling constants are given by 

R = r + (n + 2)u L- d V 1 (Alia) 

C/ = «o-(n + 8)^i- rf ^- 1 (Allb) 

^(r + q 2 +6|<7| CT ) 

These are the finite-size corrections to the bulk coupling constants r$ and uq, which are necessary for the evaluation 
for various thermodynamic quantities. 



APPENDIX B: FINITE-SIZE SCALING BEHAVIOR OF THE EVEN MOMENTS OF THE ORDER 

PARAMETER 

By definition the 2p-th moment of the order parameter of an O(n) model is given by 

< M * >»= • (B1) 

Changing the variable of integration to </? = (J7L d ) 1 / 4 <i> and by introducing the scaling variable z = RL d / 2 / ' \/Tj the 
above expression can be rewritten in the form 

<M 2p>n =(UL d ) 2 ^-^f ; 5 ; 4 . (B2) 



Using the identity p& 



rVW-Tv, . (2ffl-^r W exp (g) (-L) 



(B3) 



where D p (z) are the parabolic cylinder functions, the above expression can be rewritten in a very simple form 



r[n/2] D_„ /2 (z/V2) 
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Using now the asumptotics of D p {z) it is straightforward to obtain the asumptotic behavior of the above moments 
for i) z 3> 1 and ii) z <C 1. 
i) z ^> 1. Then one has 



<M 2p > n = {UL d /A)~ 



T[p + n/2] 



T[n/2 

ii) z <C 1. For this case the corresponding result is 

[UL d 



p(n +p+l) 



< M 2p >„ 



V 4 



£[f+f] 
r[f] 



1 2 ' 4 J 

r[f] 



r[f 



nh + j] /r[i + f] r[i + f] r[i 



r[f] 



n , p- 
2 1 4^2- 



r[i + j + §] 
2r[f + f] 



m 2r[f + i] 

+ 0(z 3 ) . 



r[f + §] 



For the susceptibility (p = 1) the above expression can be written in the following very simple form 

< Mi >„= a n + zb n + z 2 c n + 0(z 3 ), 

where a n = T(f + ±)/r(f ), 6„ = a* - % c n = a n (b n + 1/4). 

^From ( |B5| ) it follows that the asymptotic behavior of the Binder cumulant is 



1/2 
B n (z) = 1 - - 1 + - 
J \ n 



1 ^T-^l4 



wherefrom one has B n (oo) = |(1 — 1/n). 



(B5) 



(B6) 



(B7) 



(B8) 



APPENDIX C: FINITE-SIZE RESULTS FOR THE PHYSICALLY IMPORTANT CASE: d + a = 6 



In this Appendix we will report some results for the important case d + a = 6, which models the van der Waals 
type potential. Note that because of the condition d + a = 6 one now has only one independent variable, i.e. setting 
d = 4 — e directly leads to a = 2 + e. If one performs now e-expansion on the er-dependent terms one will in fact 



change the spectrum of the system from such one, where q° 



is considered as a perturbation to the short-range 



contribution (proportional to q ), to one in which q a is replaced by q + eq lng, i.e. where the long-range portion 
of the interaction will represent already a leading-order term. This is not the type of systems we are interested in. 
Therefore, in order to avoid this problem, in all the calculations below we perform e-expansion only on the d-dependcnt 
terms and retain the full e-dependence in all terms where it is stemming from the a-de pend ence of the quantities 
involved. Following this way of acting we obtain that in the case d + a — 6 the expression ( |3.21 ) for R transforms into 



RL 2 



y 1 1 + 7. ln v 



2n 

— Sa 



Scaling (yA)~bL-c \l + y 



bL- £ y{\ny -2\nL) 
d 



dy 



scaling 



(Ci) 



showing that there is an additional InL correction to the finite-size scaling theory. Definitely, keeping the terms 
proportional to L~ e one goes beyond the precision kept in the remaining part of the above equation. In accordance 
wit h th e remarks made above note that wile one does not perform an expansion of L~ e in terms of e all the terms 
in (CI) proportional to L~ e are simply corrections to scaling. But, once one performs that expansion, because of the 
lnL proportionality, these terms produce a leading-order contribution, which is quite unphysical. We believe that this 
is an artifact of the e expansion. Such a procedure (keeping the full e-dependence in some expressions) has been used 
in reference |37|] in the analysis of the scaling properties of quantum systems at low temperatures. We hope that the 
above problems can be removed by performing, e.g., a field theoretical method based on minimal renermalization at 
fixed space dimensionality |l^,^8|. This is out o f the scope of the current article. 
For the coupling constant U, instead of ( [3.2S ) one obtains 
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UL £ 



l + -(l + lny) 



?i + 8 



bL~ e SA [Lay - 2 In L] 



n + 8 



dy 



Scaling (yA)-bL- 



dy dy 2 



I 2 

scaling 



(y, 4) 



(C2) 



In this limit an additional InL correction shows up. This expression is finite in the limit y = 0. At the critical point 
it transforms into 



UL £ = 



n + 8 



-S 4 





e 2 




bSi 




n + 8 



-C-2\nL - 1 
2 



5 e 2 



2n- 



-MT e C(3). 



(C3) 



The explicit appearance of In L will affect the result of the susceptibility, which will depend upon an additional In L 
at the critical point T = T c . 

At the fixed point, for the 'characteristic' variable z, we obtain 



RL 2 



VuT e 



fixcdpoint 



71 + 8 
71 + 2 



n — 4 
n + 8 



In y 



77+2 i 4 — 77 

ebL- £ ——y + ebL- £ —-y(\ny - 21nL) 

77 + 8 2 77 + 8 



£5*4 



n - 



.,i, llg (2/>4) ///. ( 1 • !l— ) 



c9 



(y,4) 



-—eSm 
2 y 



(9 



rfcaling (y,4)-6L" £ 2— + — 



d d 2 



dy dy 2 



scaling 



(y, 4) 



(C4) 



A comparison between ( |3 ■ 28[ ) obtained for the case d + a < 6 and ( C4 ) show s that an additional In L appears in 
the expression of the variable z(y), however this does not alter the result (3.29) for z*(0), i.e z(y) evaluated at the 
critical point T — T c . In this case the term proportional to InL vanishes as we take the limit y — * 0. The result (C4) 
shows, in this way, that the Binder Cumulant at th e critical point does not depend on lni . 

Far away from criticality the susceptibility ( |3.34 ) found for the case d + a < 6 turns into 



X = Xc 



1 + e 6 ^!^ 1 (lnxoo + 5£S 4 C(3)y- 3 ) 



(C5) 



for the case d + cr = 6. Remark that the susceptibility conserves the same features as that of the case discussed in the 
body paper. 
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